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ABSTRACT 

00 . 

Context. Phoebe is one of the irregular satellites of Saturn; the images taken by Cassini-Huygens 
spacecraft allowed us to analyze its surface and the craters on it. 
Q_4 ' Aims. We study the craters on Phoebe produced by Centaur objects from the Scattered Disk (SD) 

and plutinos escaped from the 3:2 mean motion resonance with Neptune and compare our results 
with the observations by Cassini. 

Methods. We use previous simulations on trans-Neptunian Objects and a method that allows us 
2 ■ to obtain the number of craters and the cratering rate on Phoebe. 

C/3 , Results. We obtain the number of craters and the greatest crater on Phoebe produced by Centaurs 

in the present configuration of the Solar System. Moreover, we obtain a present normalized rate 
of encounters of Centaurs with Saturn of F = 7.1 x 10~" per year, from which we can infer the 
^ ' current cratering rate on Phoebe for each crater diameter. 

OO ' 

Conclusions. Our study and the comparison with the observations suggest that the main crater 
features on Phoebe are unlikely to have been produced in the present configuration of the Solar 
System and that they must have been acquired when the SD were depleted in the early Solar 
System. If this is what happened and the craters were produced when Phoebe was a satellite of 
Saturn, then it had to be captured, very early in the evolution of the Solar System. 

Key words, methods: numerical - Kuiper Belt: general 
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s-h 1. Introduction 

Phoebe is one of the irregular satellites of Saturn. It has a retrograde orbit, which suggests that it was 
captured by Saturn instead of being formed "in situ" (e.g Pollack et al. 1979). Moreover, Phoebe's 
composition is close to that derived from bodies such as Triton and Pluto, and it is different from 
that of the regular satellites of Saturn supporting Phoebe's origin as a captured body of the outer 
Solar System (Johnson & Lunine 2005). 

On 1 1 June 2004, Cassini-Huygens spacecraft encountered Phoebe in a fly-by within 2000 
km of Phoebe's surface. This encounter allowed Cassini to analyze Phoebe's surface in detail 
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which improved the previous Voyager data. Buratti et al. (2004) analyzed and characterized the 
physical properties of the surface from photometric data from Cassini VIMS (Visual and Infrared 
Spectrometer), concluding that it is rough and dusty, perhaps from a history of out-gassing or a 
violent collisional history suggested by Nesvorny et al. (2003 ). Hendrix & Hansen (120081 ) analyzed 
the first UV spectra of Phoebe with the Cassini Ultraviolet Imaging Spectrograph (UVIS) during 
the Cassini spacecraft fly-by and detected water ice. Using VIMS data, Cruikshank et al. (2004) 
concluded that Phoebe's surface is rich in organics, which is compatible with the low albedo of 
Phoebe. Porco et al (120051 ) calculated a mean density of Phoebe of 1.63gr/cm 3 from calculation 
of the volume and the determination of the mass from tracking the spacecraft. If Phoebe's sur- 
face was a mixture of rock and ice, the previous density is compatible with a porosity lower than 
~ 40% (Porco et al. 2005). Johnson and Lunine (2005) analyzed the relation between composition 
and probable porosity of Phoebe and they found that if Phoebe was derived from the same com- 
positional reservoir as Pluto and Triton, Phoebe's measure density is consistent with a porosity of 
~ 15%. 

Giese et al. (2006) presented the results of a photogrametric analysis of the high-resolution 
stereo images of Phoebe. In particular they obtained a mean figure radius of 107.2 km and a digital 
terrain model of the surface reveals significant morphological detail. The images revealed that 
Phoebe basically exhibits simple crater shapes with the only exception of the greatest impact crater 
Jason with a diameter of ~ 100 km. Several of the smaller craters present pronounced conical 
shapes which could indicate the presence of porous, low compacting material on the surface of 
Phoebe. 

Kirchoff & Schenk (2010) reexamined the impact of crater distribution of the mid-sized sat- 
urnian satellites. For Phoebe they found that the crater size frequency distribution has relatively 
constant values for crater diameters D < 1 km, but then it has a sudden and confined dip around 
D ~ 1.5 km. Beyond this dip, the crater size frequency distribution has a slow increase. This 
behaviour is unique in the saturnian satellite system and it is probably connected with Phoebe's 
origin. 

Zhanle et al. (2003 ) write a previous paper that calculates cratering rates in the satellites of the 
outer planets. They used impact rates on the giant planets obtained by Levison & Duncan ( 1997) 
and independent constraints on ecliptic comets. Their results will be compared with ours. 

As we have seen, the origin of craters on Phoebe is unclear. But the main population of objects 
that can produce craters on Phoebe are Centaurs, since they are the small body objects that cross 
the orbit of the giant planets, in particular the orbit of Saturn, and then its satellites. 

Centaurs are transient bodies between their source in the trans-Neptunian population and the 
Jupiter Family Comets. They come mainly from a sub population in the trans-Neptunian zone, the 
Scattered Disk Objects (SDOs). The SDOs are bodies with perihelion distances q greater than 30 
AU and smaller than ~ 39 AU that can cross the orbit of Neptune and eventually evolve into the 
giant planetary zone, crossing the orbits of those planets, and then the orbits of their satellites (Di 
Sisto & Brunini 2007 Levison & Duncan 1997 ). The secondary source of Centaurs are plutinos 
and the low eccentricity trans-Neptunian objects (Di sisto et al. 120101 Levison & Duncan fl 9971 ). 
Plutinos are those trans-Neptunian objects located in the 3:2 mean motion resonance with Neptune 
at a ~ 39.5 AU. They are "protected" by the 3:2 mean motion resonance with Neptune but some of 
them are long term escapers that are presently escaping from the resonance (Morbidelli 1 1 9971) . In 
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this paper we will study the production of craters on Phoebe from Centaur objects from SDOs and 
plutinos escaped from the 3:2 mean motion resonance with Neptune, as the two main populations 
of impactors. We use here previous simulations on trans-Neptunian Objects (TNOs) and a method 
that allows us to obtain directly the cratering rate on Phoebe. This study, and the comparison with 
the observations of Cassini images may help us to determine the origin of crates on Phoebe, in 
order to determine the history of Phoebe's surface and also constrain its origin. 

2. The number of SDOs 

Cratering rates depend on the number and sizes of the impactor population. Thus we must know 
the real initial number of SDOs to calculate the total number of collisions on Phoebe. Then, we are 
going to estimate the total number of present SDOs. 

f Parker & Kavelaars d20 1 Oat re-characterized the orbital sensitivity of several published pencil- 
beam surveys. They found that these surveys were sensitive to distant populations like SDOs and 
Sedna-like objects. Using this result Parker & Kavelaars (2010b) derived new upper limits on those 
distant populations and used this new limits to obtain the number of SDOs. To do this they per- 
formed a model that considered two laws for the radial distance distribution of SDOs. On the one 
hand they took a radial distance distribution of SDOs oc r~ L5 and obtained a maximum population 
of N(d > 100 km ) = 3.5 x 10 5 . On the other hand they took a uniform radial distance distribu- 
tion obtaining in this case a maximum population of N(d > 100 km ) = 25 x 10 5 . In this paper 
we take the number of SDOs greater than d = 100 km equal to N(d > 100 km ) = 3.5 x 10 5 
since the considered radial distance distribution is consistent with the one obtained by Di Sisto & 
Brunini (2007). Then the total population of SDOs with diameter greater than do will be given by 
N(d > do) = 3.5 x 10 s (do/lOO)~ s , where do must be in km and s is the index of the differential 
size distribution. There are some authors who have found a single power law size distribution for 
TNOs (Petit et al.[2000 Fraser et al. 2008 ). However, other papers suggest that the size distribution 
function (SDF) of TNOs could have a break at a diameter of ~ 60 km (Bernstein et al. 2004] Gil 
Hutton et al. 120091 Fraser & Kavelaars 120091 Fuentes & Holman|200i Fuentes et al. l2009l . The 
differential power law indexes for smaller TNOs (this is d < 60 km) found by those surveys are 
S2 = 2.8,2.4, 1.9,2.5 and 2, respectively. It seems to be enough evidence for a break in the size 
distribution of the TNO population. In particular we are going to assume that this break is also 
valid for all the dynamical classes of TNOs. 

Elliot et al. (2005) accounts for the SDF for each dynamical class in the TN region. Specifically 
for SDOs they found that the differential size distribution index is si = 4.7 for the brightest objects. 
Then taking into account the assumed break in the SDF of SDOs we consider that the power 
law SDF of SDOs breaks at d ~ 60 km to an index of between 3.5 and 2.5. We analyze those 
indexes as limit cases that give higher and lower ranges for the population of SDOs and then the 
production of craters on Phoebe. The higher value of 52 = 3.5 corresponds to a population in 
steady-state (Dohnanyi 1969) which could be the case for the smallest SDOs (Gil Hutton et al. 
2009 ). Considering all this, the number of SDOs greater than a diameter do will be given by 
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Fig. 1. Cumulative number of SDOs and Plutinos according to the size distribution laws described 
in the text. 

where C = 3.5 x 10 5 100 iI-1 (60)* 2- ' tI by continuity for d = 60 km, s, = 4.7 and s 2 = 2.5 and 3.5. 
This law is plotted in Fig. 1 with the two breaks considered. 



3. SDO collisions on Phoebe 

To study the collisions of SDOs on Phoebe and the contribution of that population to the cratering 
history of the satellite we use some of the outputs of the numerical simulation performed in a 
previous paper by Di Sisto & Brunini (2007). In that work we integrated numerically 1000 objects 
from the SD (95 real + 905 fictitious) and studied their evolution in the Centaur zone under the 
gravitational action of the Sun and the four giant planets. The computations were followed for 4.5 
Gyr, or until the test body collided with a planet, was ejected, or entered the region inside Jupiter's 
orbit (r < 5.2 AU). In that paper we also stored in a file the encounters of the fictitious SDOs with 
the planets and registered the time of the encounter, the minimum distance to the planet (q) and the 
relative velocity at this distance (v(q)). From these data we can calculate, the number of encounters 
with Saturn within the Hill's sphere of the planet. Using the particle in a box approximation and 
assuming that the geometry of the encounters is isotropic it is possible to calculate the number of 
collisions on Phoebe (N c ) through the relation: 



N c 



v. Ri 



(2) 



N e v(R) R 2 

Where N e is the number of encounters with Saturn inside its Hill's sphere of radius R, R p is the 
radius of Phoebe, v(R) is the mean relative encounter velocity of SDOs when entering the Hill's 
sphere of the planet and v, is the collision velocity of SDOs on Phoebe. v(R) can be calculated from 
v(q) registered in our outputs from the relation: 

1 1 

<■/ 

where G is the constant of gravitation and m is the mass of Saturn. 

The collision velocity on Phoebe is computed assuming that the geometry of collisions is 
isotropic and then: 

Vj = yj 



v(R) = v\q) + 2Gm(i - -) 



(3) 



ivj + v 2 Q 



(4) 
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Table 1. Radius R p , density p p , orbital velocity v p and superficial gravity times the mass gm p of 
Phoebe. Mean relative velocity of SDOs when entering the Hill's sphere of Saturn v(R), when they 
intersect the orbit of Phoebe vq and when collide with Phoebe v,-. 



R p [km] 


107.2 


P P [gr cm" 3 ] 


1.634 


v p [kms~'] 


1.71 


gm p [km s~ 2 ] 


0.5532 


v(R) [km s-'] 


4.06 


v [km s~'] 


4.65 


v, [km s~'] 


4.96 



where v p is phoebe's orbital velocity and vo is the mean relative velocity of SDOs when they cross 
the orbit of Phoebe. This velocity was computed in the same way as v(R) appropriately using Eq. 
(0 from v(q) registered in our outputs. All the mentioned velocities and the radius and orbital 
velocity of Phoebe are shown in Table 1 . 

Eq. (|2]i provides the number of collisions on Phoebe in relation to the number of encounters 
with Saturn that we had registered in our previous run. 

There are many papers based on theoretical and observational work which argue that the initial 
mass of the trans-Neptunian region was ~ 100 times greater than the present mass and decay to 
nearly its present value at most in 1 Gyr (see e.g. Morbidelli et al. l2008t . The simulation of Di Sisto 
& Brunini (2007 ) studies the evolution of SDOs in the present configuration of the Solar System; 
that is, when the SD is expected to have roughly reached its present mass and dynamic state, ~ 3.5 
Gyrs ago. Then, we can estimate the total number of collisions on Phoebe during the last ~ 3.5 
Gyrs rescaling Eq. (0 in order to account for the total SDOs population. 

Of the 1000 initial particles of our previous simulation (Di Sisto & Brunini 120071 1. 368 under- 
went 10 257 encounters within Saturn's Hill sphere. Therefore, the total number of encounters with 
Saturn of the whole SDO population in the present configuration of the Solar System is estimated 
as 

M j 368 \ / 10257 \ M 

N « = (lOOo) (-368") (5) 
where N is the number of different SDOs which have existed, in the last 3.5 Gyrs and which can 
be obtained from Eq. Here we assume that the present number of SDOs is roughly the same as 
it was 3.5 Gyrs ago. Consequently, the total number of encounters with Saturn of the whole SDO 
population through the last 3.5 Gyrs depending on the diameter is given by: 

N el (d>d ) = ^^^N(d>d ) (6) 

From this equation and Eq. (0 the total number of collisions of SDOs on Phoebe over the last 3.5 
Gyrs depending on the SDO's diameter, is given by: 

V' R^ 

N ^ d > *) = N ^ d > d o) ( ? ) 

v(R) R 

Table|2]shows some values of N c for certain diameters of the impactors. 

Depending on the values of S2, the diameter of the largest SDO impactor onto Phoebe during 
the last 3.5 Gyrs has been calculated as ranging from 110 mts to 1 .36 km. 
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Table 2. Number of collisions of SDO impactors on Phoebe with diameters d > do that produce 
craters with diameter D > D (N C (D > D a )). 



do[km] 


D [km] 


N C (D > Do) 


0.081 


1 


2-1180 


0.445 


5 


0- 16 


0.969 


10 


0-2 


3.48 


30 






4. Craters on Phoebe by SDOs 

The estimation of the size of a crater produced by a particular impact has been extensively studied. 
Schmidt & Housen ( 1987 ) present a set of power-law scaling relations for the crater volume based 
on laboratory experiments that simulate crater formation and point-source solutions. Holsapple 
(1993 ) also describe the scaling law for impact processes in a review work. The derived scaling 
laws allow us to link impacts of different sizes, velocities and superficial gravity and then obtain the 
size of a crater produced by a collision on a solar system body. Holsapple & Housen (2007 ) present 
the updated scaling laws for cratering in a recent work dedicated to interpret the observations of 
the Deep Impact event. These cratering laws are used here to calculate the craters on Phoebe. Thus, 
the diameter D of a crater produced by an impactor of diameter d can be obtained from the general 
equation (Holsapple & Housen 120071) : 
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ter- 


f Y ] 
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fe) 
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p, being the target density, g its superficial gravity, Y its strength, p, the density of the impactor 
and V; the collision velocity. This impact cratering scaling law depends on two exponents, fi and v, 
and a constant, K\, that characterize the different materials. The first term in the square brackets 
is a measure of the importance of gravity in the cratering event and the second is a measure of 
the importance of the target strength. Thus, if the first term dominates on the second term, the 
crater is under the gravity regime, and if the second term dominates we have the strength regime. 
The partition between the two size scales of impacts depends on the size of the event (Holsapple 
119931 ). Eq. (O is a convenient empirical smoothing function to span the transition between the 
gravity regime and the strength regime (Holsapple, 1993). Since Phoebe is a small satellite with a 
relatively low gravity, the strength regime can be important for the smaller craters. 

As Phoebe's density (1 .63gr/cm 3 ) is similar to sand and lower compacting material is found on 
its surface, we adopt K\ = 1.03, n = 0.41 and v = 0.4 that correspond to sand or cohesive soil in 
Holsapple & Housen (2007 ). This value of fi corresponds to materials with a porosity of ~ 30-35% 
(Holsapple & Schmidt 1987) which is compatible with the ranges of Phoebe's predicted porosity. 
The value for dry soils from Holsapple ( 1993 ), i.e. Y — 0.18 mpa, is used for the strength. 

The calculated densities of TNOs vary considerably from ~ 0.5- ~ 3 gr/cm. Although a 
dimension-density trend has been suggested (Sheppard et al. 120081 Perna et al. 120091 ). more data are 
required to confirm it. In addition; as crater experiments do not include variations in the impactor 
material, there is no data to precisely determine the dependence on the impactor density (Schmidt 
& Housen 1987 Housen & Holsapple 2003). Therefore, we assume p, = p,, which is also between 
the ranges of calculated densities in the trans-Neptunian region. By taking all this into account, 
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the diameter of a crater on Phoebe for a given impactor diameter can be calculated from Eq. © 
through: 



D = 1.03 



n -0.17 

Y 



Ptvf 



d (9) 



This equation describes simple bowl-shape craters but, as mentioned in the introduction, Cassini 
images of Phoebe reveal basically simple crater shapes with the only exception of Jason with a 
diameter of ~ 100 km (Giese et al. 2006 ). Hence we use Eq. (01 for calculating the diameters of all 
craters on Phoebe without any further correction for transient-to-final size. By combining Eq. (Q 
and (|9]l it is possible to calculate the number of craters on Phoebe according to the diameter of the 
crater. Fig. ((2j shows the cumulative number of craters, with diameters greater than a given value 
for the two size distribution power laws for smaller SDOs on Phoebe. Note that the different slopes 
in the number of craters for each curve is due to the difference in both indexes S2 considered. As 
we mentioned before, there is a limit impactor diameter that accounts for the transition between 
the gravity regime and the strength regime. This diameter can be obtained equating the first and 
second terms of Eq. (0. This limit impactor diameter is di = 2>(flmts, which produce a limit crater 
of Di = 4.2 km. Thus, for crater diameters D < D\, the production of craters on Phoebe is under the 
strength regime and for D > D\, the production of craters is under the gravity regime. Note that D\ 
depends strongly on the assumed value of strength which is actually unknown. Then considering 
other values of the strength for less cohesive soils as terrestrial dry desert alluvium of Y — 65 kpa 
(Holsapple & Housen 120071 1 and Surface Lunar regolith Y = 10 kpa (Holsapple l2011l i. D\ could 
take the values 1.5 km and 233 mts respectively. We will consider in the following that D\ = 4.2 
km but it must be taken into account that D\ can be as small as 233 mts. 

Since in the strength regime the crater diameter depends linearly on the impactor diameter, the 
relation between the cumulative number of craters on Phoebe and the crater diameter follows the 
same power law relation as that followed by the number of SDOs. For D < 4.2 km the cumulative 
number of craters on Phoebe follows a power law with a cumulative index of 1 .5 and 2.5, according 
to the value of 52 = 2.5 or S2 =3.5 respectively. For D > 4.2 km, this is in the gravity regime, 
the crater diameter does not depend linearly on the impactor diameter. Therefore, we fit a power 
law for the cumulative number of craters on Phoebe depending on the crater diameter of index 
2.8. Kirchoff & Schenk (2010) found that the crater size frequency distribution for Phoebe has a 
cumulative index of 2.348 for D - 0.15 - 1 km and 1 for D — 1 - 4 km. We can see that, for very 
small craters, this index is very similar to our value of 52 = 3.5 or cumulative index of 2.5. This 
is consistent with the fact that the size distribution of very small objects is expected to approach a 
Donhanyi size distribution (S2 = 3.5) and then the craters produced by those small projectiles (that 
are in the strength regime) have to follow the same power law size distribution. Besides Kirchoff 
& Schenk (2010) found that for D = 1 - 4 km, Phoebe's crater distribution has a shallow slope and 
this implies that Phoebe has a deficiency of craters with D ~ 1 .5 km. This change of slope cannot 
be explained by our method and our proposed contribution of Centaurs from the SD to the craters 
on Phoebe, unless that the considered SDF of SDOs is different. Anyway more work in relation to 
another source of craters on Phoebe, as planetocentric objects, is needed and also its connection to 
the origin of the irregular satellite itself. 

According to the differential size distribution index S2, the largest crater on Phoebe produced 
by a Centaur from the SD has a diameter of between 1.4 km and 13.5 km. Table [2] shows the 



8 



R. P. Di Sisto & A. Brunini: Origin of the craters on Phoebe 



o 




0.01 0.1 1 10 

D [km] 

Fig. 2. Cumulative number of craters with diameters greater than Do in the last 3.5 Gyrs, produced 
by Centaurs from SDOs on Phoebe. Filled line corresponds to the differential power law index 
S2 — 3.5 and the dotted line to 52 = 2.5 

cumulative number of craters on Phoebe greater than certain diameters produced by Centaurs from 
the SD in the current configuration of the Solar System during the last 3.5 Gyrs. Since the largest 
crater on Phoebe has a diameter of ~ 100 km, it is unlikely that it was produced by a recent collision 
of an SDO. This will be discussed in a following section. 

5. Rate of SDO collisions on Phoebe 

From our outputs we can calculate the number of encounters within the Hill's sphere of Saturn as a 
function of time. In Fig.[3]is plotted the normalize cumulative number of encounters as a function 
of time. The whole plot can be fitted by a log-function given by /(f) - a +b log t, where a = -3.24 
and b = 0.19. The total cumulative number of encounters with Saturn for all diameters for each 
time can be obtained from the plot and Eq. ||6). We have calculated the number of collisions on 
Phoebe from the number of encounters with Saturn from our outputs (Eq. [7]), and this number 
can be calculated for each time in the integration. Then the cumulative number of collisions (or 
cumulative number of craters) on Phoebe for all diameters as a function of time can be obtained by 
multiplying the fraction of encounters obtained from Fig.[3]by the number of collisions N c (d > do). 

As we can see from Fig.[3]the rate of encounters and then the rate of collisions on Phoebe was 
high at the beginning but it has been decreasing up to the present. The first Myrs the shape of the 
curve is purely arbitrary due to initial conditions but then it will be significant and is stabilizing. In 
the last ~ 3.5 Gyrs the rate has been almost constant. In fact, it is possible to fit a linear relation 
to the last 3.5 Gyrs of Fig. [3]given by g(t) = Ft + c, where F = 7.1 X 10~ n and c = 0.69. This 
linear approximation allows us to calculate a present rate of craters for a given diameter. The slope 
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Fig. 3. Fraction of encounters of SDOs with Saturn. The linear fit to the data the last 3.5 Gyrs is 
also shown. 

of this linear function F = 7.1 x 10" 11 is the present normalized rate of encounters of SDOs with 
Saturn per year. To obtain the present rate of encounters with Saturn for each diameter we must 
multiply F by N er (d > do) (obtained from Eq (O). And then the present rate of collisions on Phoebe 
for each diameter can be obtained multiplying F by N c (d > do) (obtained from Eq (0). Similarly, 
the current rate of cratering on Phoebe for craters greater than a given diameter can be obtained 
from the current rate of collisions and the relation ^ between the diameter of the impactor and 
the diameter of the crater. Thus, for example the current cratering rate on Phoebe of Centaurs from 
SDOs that produces craters with D > 1 km is between 1.4 x 1CT 10 and 8.3 x 10~ 8 craters per year 
(depending on the 52 value); this is at least ~ 80 craters with D > 1 km in the last Gy. The current 
cratering rate of craters with D > 5 km is at least 1.14 x 10~ 9 craters per year, this is ~ 1 crater 
with D > 5 km in the last Gy. 

Zhanle et al. ([2003) calculate cratering rates in the satellites of the outer planets. They obtained 
a cratering rate on Phoebe for craters with D > 10 km of 8.6 x 10 year . In this study we obtain 
a cratering rate on Phoebe of craters with D > 10 km of between 2.7 x 10~ 12 and 1.4 x 10~ 10 craters 
per year, depending on S2- Our value for S2 = 3.5 is very similar to the previous calculation of 
Zhanle et al. (120031 

6. The contribution of escaped Plutinos to the craters on Phoebe 

Plutinos could be another source of craters on Phoebe. In a recent paper, Di Sisto et al. (2010) 
study the post escape evolution of Plutinos when they escaped from the 3:2 mean motion resonance 
with Neptune, and in particular their contribution to the population of Centaurs. In that work, we 
perform two sets of numerical simulations in order first to identify the plutinos that have recently 
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escaped from the resonance and second to follow their evolution under the influence of the Sun and 
the four giant planets. This numerical simulation considers the evolution of plutinos in the present 
configuration of the Solar System, this is as we did for SDOs, when the trans-Neptunian region is 
expected to have reached roughly its present mass and dynamical state, ~ 3.5 Gyrs ago. Following 
the same analysis that we made for SDOs described in Sects. [3]and[4] we calculate the number of 
craters on Phoebe from escaped plutinos the last 3.5 Gyrs, and also the greater impactor and crater. 
In the numerical simulations by Di Sisto et al. ( 120 1 Oi l, we start with 20 000 initial massless particles 
in the 3:2 mean motion resonance, 671 of them undergone 20459 encounters within Saturn's Hill 
sphere during the integration. We obtained that the mean relative encounter velocity of plutinos 
when entering the Hill's sphere of Saturn is v(R) = 4.57 km/s, the mean relative velocity of plutinos 
when they intersect the orbit of Phoebe is vo = 5.12 km/s and the collision velocity of plutinos on 
Phoebe is v, = 5.4 km/s. We take the present number of plutinos from de Elfa et al. (20083, but 
consider that the size distribution breaks at d ~ 60 km with the two limit power law indexes S2 = 3.5 
and 2.5 as we adopted for SDOs to be consistent (see Sect. |2j. Those indexes give the highest and 
lowest number of SDOs and then the highest and lowest production of craters on Phoebe. Then the 
present cumulative number of plutinos is given by: 



where C = 7.9 x 10 9 (60)''~ 3 by continuity for D = 60 km and the cumulative power law index p 
adopt the values 2.5 and 1.5 (p = s - 1). This law is plotted in Fig 1. 

Considering all this and our method described above, we obtained that the greater escaped 
plutino impactor on Phoebe has a diameter between 1.5 mts and 102 mts that produces a crater 
between 19.3 mts and 1.3 km respectively, depending on the power index p of the size distribution 
of plutinos. Also, we can obtain the number of craters on Phoebe from escaped plutinos. We have 
at least two craters greater than 1 km on Phoebe from plutinos. Comparing this with the values 
obtained for the contribution of SDOs, it can be stated that the number of craters produced by 
escaped plutinos on Phoebe is negligible with respect to the contribution of SDOs. Also it can be 
stated that the greater craters are signed by the contribution of SDOs. 

7. Discussion 

In the previous sections, we have calculated the production of craters on Phoebe considering the 
present population in the SD and Plutinos. However - as mentioned - there are many papers based 
on theoretical and observational work that argue that the initial mass of the TN region was ~ 100 
times greater than the present mass (see e.g. Morbidelli et al. 2008 ). Observations predict a current 
mass of the Kuiper Belt that is very small with respect to that required for models to grow the 
objects that we see. The mass depletion due to a strong dynamical excitation of the Kuiper Belt 
is thought to be the scenario for this "mass deficit problem". There were several models that try 
to describe the mass depletion; the last model that described this mechanism is the "Nice model" 
where the Kuiper Belt had to be significantly depleted before the time of the LHB (Levison et 
al. 2008 ). The "Nice Model" assumes the giant planets initially in a more compacted region from 
~ 5.5 to ~ 14 AU and a planetesimal disk of a total mass of ~ 35 Mr that extends beyond the orbits 




for 



for 



D > 60km, 



D < 60km, 



(10) 
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of the giant planets up to ~ 34 AU. The interaction between the planets and planetesimals makes 
the giant planets migrate for a long time removing particles from the system. After a time ranging 
from 350 My to 1.1 Gy, Jupiter and Saturn cross their mutual 1:2 mean motion resonance. Then, 
the eccentricities of Uranus and Neptune drives up and those planets penetrate into the planetesimal 
disk. This destabilizes the full disk and the planetesimals are scattered all over the Solar System. 

Beyond the model and the mechanism responsible for the mass depletion of the trans-Neptunian 
zone, we can consider that primitive SDOs (that were 100 times more numerous than the present 
ones) follow the same dynamical evolution than the present population when they enter the plan- 
etary zone as Centaurs. Then we can calculate in the same way as we did in the previous sections 
and with the same model the cratering on Phoebe assuming an initial population of SDOs 100 times 
the present population. 

This is an estimation since we have to know the real initial scenario of formation of the Solar 
System and in particular of SDOs. However when a SDO enter the Centaur zone, inside the orbit 
of Neptune, its dynamical evolution is governed by the giant planets and then the particular initial 
scenario can be considered secondary for the present study. 

Doing that, we obtain that the greater impactor on Phoebe during the age of the Solar System 
has a diameter between 2.4 to 8.6 km and produces a crater of 21.6 to 64.2 km. The value corre- 
sponding to S2 = 3.5 (64.2 km) is in agreement (within the expected errors and statistical fluctua- 
tions) with the observation of the great crater "Jason" on Phoebe with a diameter of ~ 100 km. The 
number of craters greater than a given diameter can be obtained increasing 100 times the values 
obtained in Sect. [4] 

Recently, Cassini images of Phoebe have allowed to study its surface and the craters in it. 
Kirchoff & Schenk (2010) obtain, a cumulative crater density for D > 5 km of 2233 + 1117 based 
on crater counting from Cassini images. From our model, assuming an initial number of SDOs 100 
times the present population, we obtain N C (D > 5km) = 12 - 1640, again in a good agreement with 
the values obtained by Kirchoff & Schenk (120101) S2 = 3.5. 

8. Conclusion 

We have studied the production of craters on Phoebe from SDOs and escaped plutinos that have 
reach the Saturn zone in the present configuration of the Solar System. We have obtained that the 
contribution of escaped plutinos is negligible with respect to the contribution of SDOs. We have 
obtained that both the number of craters and the greater crater on Phoebe produced by SDOs cannot 
account to the observations. But if we take into account that the initial mass of the trans-Neptunian 
region was 100 times the present one, we match the craters produced by SDOs on Phoebe with the 
observed characteristics of the satellite if S2 = 3.5. 

Those considerations imply that the main cratering features of Phoebe must be acquired when 
the SD had being depleted at the early times of evolution of the Solar System. More than that, if 
the "Nice model" describes correctly the scenario of the origin of the Solar System, the scatter- 
ing inward of planetesimals by Neptune and Uranus in that model must be similar to the present 
scattering in our model, and the TNOs have to lose memory when they arrive to Saturn. 

If this is what happened and the main crater characteristics on Phoebe were produced when 
Phoebe was a satellite of Saturn, the great agreement of our model with the observations constrain 
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the time when Phoebe had to be captured, very early in the evolution of our Solar System. This was 
also suggested by Levison et al. (2008 ). 

We have obtained that the present normalized rate of encounters of SDOs with Saturn is: F = 
7.1 x 10" 11 per year. From this number we could obtain the present cratering rate on Phoebe for 
each crater diameter. 

We have compared the size crater distribution on Phoebe obtained from our model with the 
observations of craters by Kirchoff & Schenk (2010). Our crater size frequency distribution agree 
with that obtained by Kirchoff & Schenk for very small impactors that produce craters with D = 
0.15 - 1 km. This distribution follow a power law with a cumulative index of 2.5 consistent with a 
Donhanyi size distribution. For craters of D — 1-4 km Kirchoff & Schenk (2010) found a shallow 
slope and a deficiency of craters with D ~ 1 .5 km. This change of slope cannot be explained by our 
method and the contribution of Centaurs from the SD. More work in relation to another source of 
craters on Phoebe, as planetocentric objects is needed, and also its connection to the origin of the 
irregular satellite itself. 
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